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We consider the Langevin dynamics of a semi-dilute system of chains which are random polyam- 
pholytes of average monomer charge q and with a fluctuations in this charge of the size Q" 1 and 
with freely floating counter-ions in the surrounding. We cast the dynamics into the functional 
integral formalism and average over the quenched charge distribution in order to compute the dy- 
namic structure factor and the effective collective potential matrix. The results are given for small 
charge fluctuations. In the limit of finite q we then find that the scattering approaches the limit of 
polyelectrolyte solutions. 



I. INTRODUCTION 
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The properties of polyamphoJ^jtes are exceedingly complex and have recently been the subject of various theoretical 
and simulational investigationsETQ ranging from single-chain properties, the properties of solutions, to the absorption 
on walls. The phase behavior of these systems is rather rich and must be described by numerous parameters!! When 
considering random, quenched polyampholytes in statistical mechanics it is well known that the disorder requires the 
implementation of of methods such as the replica tricktL 

In this paper we intend to draw upon previous worked for polyelectrolytes to describe these systems whera there is 
random and quenched charge. It is notoriously difficult to deal with such systems in statistical mechanicstj, by the 
necessity to introduce replicas. However, as has been done for spin glassesliil it is^jeey convenient to formulate the 
dynamics in the formalism of Martin, Siggia, and Rose and of Dominicis and PelititatJ, to gain results. Here we only 
restrict ourselves to the case where the dynatnics of the random polyampholyte chains are Rouse-like. However, the 
system here differs from those of spin glasses&til in that not the interaction but the charges determining its sign and 
, magnitude are the sources of randomness. An average over this quenched disorder as a rule then produces a more 

■ complex average for the disorder than in the cases mentioned above. Here we introduce such approximations, and 
] finally also compare the results to those for the dense polyelectrolyte solutions which has been addressed in a similar 

■ formalism previously!!!! 

We introduce the model and its approximations in the next section using the starting point of chains with charges 
uncorrelated along the backbones or between chains. We also retain the dynamics of pointlike ions in the fluid in which 
the solution of polyampholyte chains is to be found. The dynamics are then transformed into a form with collective 
variables after averaging over the disorder (quenched, random charge distribution), with the use of the random phase 
approximation (RPA), and by considering a specific limits of the distribution of the random charges. The structure 
factor and effective interactions are computed in the subsequent section and the results compared with those of the 
polyelectrolytes. The contribution of the counter-ions to the effective potential is the familiar Debye-Hiickel form, 
originating in usage of a quadratic collective approximation. A further contribution comes from the randomness. 
Since the RPA is being implemented (Rouse modes are being used) we are not considering here any collapsed or 
"microgel-like" phases, but work in the limit that the system is sufficiently dense to warrant the RPA. 
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The model consists of a number, N p , of chains of equal length, L, which behave according to the high friction regime 
of the Brownian dynamics. The chains are modeled by the Edwards Hamiltonian with Kuhn length b. Associated 
with each segment (each arc position s € [0, L\) of the chains is a quenched charge of magnitude qi(s)ds for the ith 
random polyion chain, which is assumed to have a Gaussian distribution, about an average charge density q. Within 
this system there is also a set of N c counter-ions, such that on average the random system is electrically neutral. 
Assuming that the counter-ions have a valence of one this is expressed as, 

= qN p L + N c . (2.1) 

Describing the random walks by variables rj(s,t) and the counter-ions, which are assumed to be structureless, by 
Xi(t), the total energy of the system is given by the Hamiltonian, 



II. THE MODEL 
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The strength of the electrostatic interaction is given in terms of the Bjerrum length, 

where e represents the electrical charge, e the dielectric constant, and k B T = 1/(3 the temperature multiplied by 
the Boltzmann constant. The polymer chains above can also subjected to an excluded volume interaction of a 
strength given by v. The distribution of the charges (?i(s) is assumed to be Gaussian for every monomer on the chain 
independently: 

Np L nr 
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Here Q indicates the magnitude of the fluctuations and Lq represents an average excess charge per macro-ion. The 
product over arc-length of the chains above is to be understood in a discretized sense. Note that this means that even 
if the parameter q were considered to be zero, the chains would only be statistically neutral. 



A. Langevin Equations 



To start with, we follow the ideas of Fredrickson and Helfand in investigating the dynamics of neutral polymerslll. 
The dynamics of the monomer positions of the two species of the system are now given by the appropriate Langevin 
equations; for the random polyions these are: 



The Langevin equations for the counter-ions are given by, 
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The constants ( p and £ c are the friction coefficients of the polyion segments and the counter-ions, respectively. 
The expressions fi lT -(s, t) and fj, x (i) represent the random, Gaussian noise exerted upon the individual particles. 
The averages of these forces are zero and the correlations are given by the familiar expressions (fi tX (t)fj jX (t')) = 
2k B T(- 1 5 ij 6(t - t')l and (f;, r (s, t)f j<r (s', *')) = 2k B T(- 1 S i jS(t - t')S(s - s')l. There are N p and iV c such equations, 
respectively. For the sake of simplicity it is helpful to rescale the lengths above using a scale of temperature k B T = 1 
and units of length b = r-. 

Now it is straightforward to implement the functional formalism of Martin, Siggia, and RosetHl (MSR), and of 
Dominicis and PelitiEj, in which the definition of the functional delta function is used to convert the differential 
equations into a functional form: 
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where J is the Jacobian of the transformation. Similarly, a generating functional can be defined schematically by 
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such that averages can be computed in terms p£ the functional derivatives with respect to h of Z. However, upon 
making use of causality in taking the averagestj the Jacobian can be taken as being equal to 1 . This property of 
the MSR-formalism is especially useful when quenched disorder is present in the problem and it has been applied 
previously to numerous problems, including the study of spin glassesEU and of manifolds in a random potentialOEZI. 
The rules based on causality entail that any averages of products containing hatted fields (i.e. the field conjugate to 
the Langevin expression) vanish if the latest time is found in the argument of the hatted field. If two times are equal, 
that of the hatted field is made infinitesimally later to produce another vanishing average. 
It is convenient to to introduce the following collective variables: 
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The collective variables with the superscript (1) are the usual collective densities of the polymer and the counterions in 
Fourier space. The collective variables with the superscripts (2) refer to conjugate dynamical variables of the density. 
In the first two instances above these are simply density expressions for the monomers. The remaining definitions are 
those of collective variables for the components of the system. The MSR functional (2.7) can now be rewritten in the 
form, 
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The sum over k is taken over all of Fourier space. Here we have introduced the matrices with the two density variables, 



(2) 



and the respective interaction matrices, 
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for the electrostatic potential, and for the excluded volume interaction: 
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The appearance of the factor i in the expressions concerning the potentials above is purely a formality and comes 
from the fact that the dynamical formalism is being used here. 

It is now straightforward to impleaaent the quenched charge average of equation (2.4). In order to facilitate this we 
introduce a conjugate vector field <^li3 such that we can write the electrostatic part of the equation as 
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The normalization of the Gaussian integration above is given by the factor Af<p. Consequently, we a re now in a position 
to take an average with respect to the quenched disorder of the expression above using eq. (2.4). The fact that the 
MSR-formalism has a normalization equal to one permits the evaluation of this quenched avera ge w ithout the need 
to resort to the replica trick. Denoting the g-averaged quantities [having used the distribution ( |2 .4j )] by an overbar 
the following result is imminent: 
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We see that the effect of disorder is to alter the interaction between the components of the system, by a term 
dependent on the configuration of the system itself. This is familiar from averages over Gaussian distributions of 
disorder. Integrating over the fields (j> leads to the disorder-averaged expression of the MSR generating functional 
such that the electrostatic part of L\ becomes replaced by, 
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Here the trace (Tr) is taken over the (k, k')-space, over the times and the density matrix indices. The definition of 
the matrix W is, 
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It is always associated with a factor of 1/Q and it can be seen that this matrix is proportional to the mass density of 
polyampholyte chains in the system. This expression (2.18) together with the remaining terms of the system given 
earlier, and which are still contained in Cq is exact, but must be approximated in order to make further progress. 



B. Approximation for small 



and finite q 



In order to be able to approximate further the expression above and to compute the dynamical structure factor of 
the present system of polyampholytes we shall assume the flu ctuat ions of the quenched charges to be small, such that 
we can expand the determinant and invert the matrix of eq. ( 2.18 ) to first order in Q^ 1 . 

Dealing first with the determinant we see that 
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The above matrix, M, is defined by: 
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The result follows from the isotropic k-integration (note the definitions of the collective variables) and the off-diagonal 
property of the interaction matrix. We have used the fact that the trace also entails taking the summations with 
5{k-k') andtf(i-i'). 

The matrix (V~ 
an expansion: 



1 + WQ 1 ) 1 sandwiched between the two collective densities in equation ( 2.1£ ) is treated also by 
WQ-T 1 = V(l + Q- 11 
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The second term above is analogous to the term sometimes denoted by "polyampholyte attraction"!!] 
definitions of the matrices involved the second term of the expansion above looks as follows: 
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This matrix is non-diagonal in the time and Fourier component variables. However, it is totally diagonal in the 
monomer labels, which follows from the uncorrelated nature of the distribution of charges. Introducing a correlation 
would be a straightforward generalization of our methods. 

We wish to implement the random phase approximation in the collective variables and deal with contributions in 
the exponent up to quadratic order in the MSR f uncti onal. Consequently, one o f the simplest approximations entails 



the substitution of the average of the expression ( 2.24 ) above into the equation ( 2.18 ) to create an effective potential 
between the collective variables. This approximation is best when the polyampholyte chains have small fluctuations 
and retain a more-or-less Gaussian distribution in the dense, yet unentangled, system considered here. For these 
purposes we consider polyions of which the excess charges are not too small to-describe the collapsed regime. The 
collapsed regime must be descibed in a completely different formalism (see, e.g.Ej). 
Let the following matrix be defined for an arbitrary j 6 {!,..., N p }: 
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We shall use th is matrix, where averaging has taken place using the distribution given with Co in approximation [see 
equation ( [2.26j )]. 

The computation is continued by the transformation to a random phase approximation (RPA) for the collective 
variables introduced earlier. By denoting the appropriate average of a functional Q by 

<Q) = N J | jf VvjVrj I j JJ VxjVki | Qe-^ , (2.26) 

such that the integration for the structure factor and response functions takes place over the collective variables p 
with mean (p) and distribution exp — ^p ■ Sq 1 ■ p, where 

S = (p(k,w)p(k',u}')) . (2.27) 

The derivations of these quantities are straightforward and follow exactly along the lines of the paper of Fredrickson 
and HelfancO. The results are summarised in the Appendix. Furthermore, at very large length-scales it can be 
assumed that the elements of the matrix w are proportional the related structure factor. 

III. DYNAMIC STRUCTURE FACTOR 

From the work in the previous section we are now in a position to write down the expression for dynamic collective 
averaging. This is given by: 

<Q} co « =Mj V P_PP_ V ^ exp (-P C (S j + A)p c - p p (S -j, + E + q 2 A)p p - pjApJ . (3.1) 
Here Af is the appropriate normalizing factor. The matrix A has been defined through the approximations above as, 

a = v +^V^ w - (3 - 2) 

Integration over the counter-ion collective coordinates now gives the result, 

(Q) coU =tfj 2>£ p Qexp (-£ p (S "j + v + g 2 A tfA(8£ + A^A)^ 

= N J VpQexp (-^(Soj, + B)pJ , (3.3) 

with B = v + q 2 A — (7 2 A(S 0c 1 + A) _1 A. The (l,l)-component of the inverse of the expression above is the dynamical 
structure factor. By denoting the interaction term above as the matrix B the structure factor G has the form: 

q _ 'S'op.ll ~ -^22'5'op,12'5'op,21 ^ q 

(1 + B2lSo p ,l2){l + Bi2So p ,2l) 

We point out that it is less tedious to compute the structure factor by taking the RPA average over the counter-ions 
before the integration over </>. 
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A. Effective Interaction 



The above form can b e ma de to look even far simpler by noticing that it is convenient to perform the average over 
the counter-ions in eq. ( 2.16| ) already. The potential expression Vk then becomes replaced by which contains 



the screening, and the matrix B above becomes simply, 

B = v + q 2 V cS -q 2 Q- 1 V ctf -W -V cS . (3.5) 

It represents the effective interactions between the poly-ionic collective variables after averaging over the disorder. 
The term V eff represents the effective interchain and intrachain interaction Coulomb interaction before the averaging 
over the quenched charges on the polyampholyte chains, after integration over the free counter-ion (collective) degrees 
of freedom. We note the typical signature of the random quenched polyampholyte interaction!!] in the appearance of 
a term oc fc -4 and proportional to Q~ 



For R g k <C 1 it is simple to express the matrices [where w is related to W as in eq. ( 2.25 )1 
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The effective Coulomb interaction as above is the same is found in the case of polyelectrolytes0 with the charge 
density along the chain replaced by q. It is easy to see that the off-diagonal terms above in the limit lo — ► return 
the familiar Debye-Hiickel screening potential. The formalism of coupled, dynamical equations leads to the additional 
diagonal component in the matrix above, which would not have resulted had the counter-ion dynamics been neglected 
and a Debye-Hiickel potential used from the start. The definition of the Debye screening parameter is given by 

k 2 = \ B Pc (3.8) 

as usual. 
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B. Structure Factor 



By substitution of all the above into the expression for the structure factor, eq. (3.4), we are now in a position 
to discuss its properties. We emphasize again that in the approximations used, terms up to order l/Q can only be 
considered, where the effect of the term Q~ 1 q 2 V c ffWV c ff oc 0{q 2 X 2 B / (Qk 2 )). We plot the typical structure factor in 
Figure 1. A peak is noticed, as for polyelectrolytes and in in the limit as Q — > oo and the results for polyelectrolytes 
are and must be recovered. For extremely smalLscattering vector k and frequency lo a behavior of the structure factor 
proportional to k 2 /lo 2 is expected as shown mli-l 

The general properties of the peak depend on the parameters q,n,Xs,v, the polyampholyte density. Firstly, it is 
noticed that the average charge fraction q per monomer when increased causes a decrease in the peak height. This is 
plotted in Figure 2. Figure 3 shows that an increase in the screening parameter, or a corresponding increase of salt 



(see eq. (3.8)) concentration of the system causes an increase of the peak, whereas a decrease of the maximal height 
is the result of turning on the excluded volume potential for .the .chain (Figure 4). These properties are in general 
agreement with those already investigated for polyelectrolytesli3EJ. 

Secondly, the role of the quenched fluctuations of polyampholyte charge, which are contain ed w ithin the parameter 
Q auch that (q 2 X 2 b /(Qk 2 )) is small need to be considered. By comparing with equation (|3.5| ) it is clear that the 
fluctuations-related term is proportional the Bjerrum length squared in distinction to the other term of the interaction. 
Consequently in Figure 5 we give a plot of three curves of the dynamical structure fact at finite frequency lo, for a 
value of Q = 10. It is seen that an increase of the Bjerrum length leads to a decrease of the maximal height. In order 
to understand the magnitude of the role of the fluctuations we now compare the polyampholyte (with Q = 10) to a 
system which is almost completely a polyelectrolyte (Q = 10 6 ) in Figure 6. A charge density of q = 0.01 has been 
assumed in these plots of the maximum at frequency lo = 0.05 in dependence of Ab, with k = 1. The higher peak 
of the polyampholyte in distinction to the polyelectrolyte (of the same backbone charge density) is clearly shown, 
although the changes in the regime where our approximations are most valid are small. 

We note that if the Bjerrum length were to be increased further, the criterion for counterion condensation X B q — 1 
would be achieved and our gaussian model would lose its applicability. 



7 



IV. DISCUSSION 



We have briefly shown how the formalism of Fredrickson and Helfand can be used to describe the collective dynamics 
of randomly charged chains, which are still assumed to be in the Gaussian regime. The structure factor shows a marked 
peak at small frequencies and low wavevector dependent on the charge fluctuation strength Q this peak can be 
understood to emerge from the mainly polyelectrolyte behavior of our system. The effects of the small polyampholyte 
quenched charge fluctuations are in increasing the height of the dynamical polyelectrolyte peak. Here the dependence 
upon the Bjerrum length are particularly evident. 

Clearly the approximations we have made to the very lowest orders are amenable to improvement, and it is true 
that a self-consistent treatment of the problem of polyampholytes will provide a better picture. This is especially so 
in respect to the Gaussian chain assumption which has gone into the computation of our collective (non-int£|Cacting) 
distribution. Since polyampholyte conformations clearly occur in many non-gaussian (e.g. collapsed) phasesEroO, the 
altered single chain properties are the future topic of our investigations. 
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APPENDIX: RPA RESULTS 

The RPA has been computed before by Fredrickson and Helfand0 and we give the results for the conditions 
kR g < 1: 

(p«(-k,-^W(k, W )) 
(pW(-k,-c)k./,( 2 )(k, w )) 
( (0 W(-k,-c)k. /9 ( 2 )(k,^)) 
All remaining pair wise correlations are zero. 
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FIG. 1. Plot of the dynamical structure factor of the polyampholyte solution case. Here Q = 5, the friction coefficients equal 
1, the Bjerrum length is set to 1 and the chain length 100, the average charge fraction 0.2, the Debye length equals 1, and there 
is no excluded volume potential. The vertical units are arbitrary. 
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FIG. 2. Plot of the dependence of the structure factor for ui = 0.05, v = 0, Ab = 1, Q = 5, k = 1, and L — 100 for three 
different values of the average monomer charge q. 




2 4 6 8 

k 



FIG. 3. Plot of the dependence of the structure factor for uo — 0.05, v = 0, Ab = 1, Q = 5, q = 0.2, and L = 100 for three 
different values of the parameter k. 
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FIG. 4. Plot of the dependence of the structure factor for ui = 0.05, q = 0.2, Ab = 1, Q = 5, k = 1, and L = 100 for three 
different values of the excluded volume interaction strength. 
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